Tracer-based lipidomics enables the discovery of disease-specific candidate biomarkers in mitochondrial ß-oxidation disorders.

Carnitine derivatives of disease-specific acyl-CoAs are the diagnostic hallmark for long-chain fatty acid ß-oxidation disorders (lcFAOD), including carnitine shuttle deficiencies, very-long-chain acyl-CoA dehydrogenase deficiency (VLCADD), long-chain 3-hydroxyacyl-CoA dehydrogenase deficiency (LCHADD) and mitochondrial trifunctional protein deficiency (MPTD). The exact consequence of accumulating lcFAO-intermediates and their influence on cellular lipid homeostasis is, however, still unknown. To investigate the fate and cellular effects of the accumulating lcFAO-intermediates and to explore the presence of disease-specific markers, we used tracer-based lipidomics with deuterium-labeled oleic acid (D9-C18:1) in lcFAOD patient-derived fibroblasts. In line with previous studies, we observed a trend towards neutral lipid accumulation in lcFAOD. In addition, we detected a direct connection between the chain length and patterns of (un)saturation of accumulating acylcarnitines and the various enzyme deficiencies. Our results also identified two disease-specific candidate biomarkers. Lysophosphatidylcholine(14:1) (LPC(14:1)) was specifically increased in severe VLCADD compared to mild VLCADD and control samples. This was confirmed in plasma samples showing an inverse correlation with enzyme activity, which was better than the classic diagnostic marker C14:1-carnitine. The second candidate biomarker was an unknown lipid class, which we identified as S-(3-hydroxyacyl)cysteamines. We hypothesized that these were degradation products of the CoA moiety of accumulating 3-hydroxyacyl-CoAs. S-(3-hydroxyacyl)cysteamines were significantly increased in LCHADD compared to controls and other lcFAOD, including MTPD. Our findings suggest extensive alternative lipid metabolism in lcFAOD and confirm that lcFAOD accumulate neutral lipid species. In addition, we present two disease-specific candidate biomarkers for VLCADD and LCHADD, that may have significant relevance for disease diagnosis, prognosis, and monitoring.

Discovery of novel diagnostic biomarkers for Sjögren-Larsson syndrome by untargeted lipidomics.

AIM: Sjögren-Larsson syndrome (SLS) is a rare neurometabolic disorder that mainly affects brain, eye and skin and is caused by deficiency of fatty aldehyde dehydrogenase. Our recent finding of a profoundly disturbed brain tissue lipidome in SLS prompted us to search for similar biomarkers in plasma as no functional test in blood is available for SLS. METHODS AND RESULTS: We performed plasma lipidomics and used a newly developed bioinformatics tool to mine the untargeted part of the SLS plasma and brain lipidome to search for SLS biomarkers. Plasma lipidomics showed disturbed ether lipid metabolism in known lipid classes. Untargeted lipidomics of both plasma and brain (white and grey matter) uncovered two new endogenous lipid classes highly elevated in SLS. The first biomarker group were alkylphosphocholines/ethanolamines containing different lengths of alkyl-chains where some alkylphosphocholines were > 600-fold elevated in SLS plasma. The second group of biomarkers were a set of 5 features of unknown structure. Fragmentation studies suggested that they contain ubiquinol and phosphocholine and one feature was also found as a glucuronide conjugate in plasma. The plasma features were highly distinctive for SLS with levels >100-1000-fold the level in controls, if present at all. We speculate on the origin of the alkylphosphocholines/ethanolamines and the nature of the ubiquinol-containing metabolites. CONCLUSIONS: The metabolites identified in this study represent novel endogenous lipid classes thus far unknown in humans. They represent the first plasma metabolite SLS-biomarkers and may also yield more insight into SLS pathophysiology.

Hepatic triglyceride content is intricately associated with numerous metabolites and biochemical pathways.

BACKGROUND AND AIMS: Non-alcoholic fatty liver disease (NAFLD) is characterized by the pathological accumulation of triglycerides in hepatocytes and is associated with insulin resistance, atherogenic dyslipidaemia and cardiometabolic diseases. Thus far, the extent of metabolic dysregulation associated with hepatic triglyceride accumulation has not been fully addressed. In this study, we aimed to identify metabolites associated with hepatic triglyceride content (HTGC) and map these associations using network analysis. METHODS: To gain insight in the spectrum of metabolites associated with hepatic triglyceride accumulation, we performed a comprehensive plasma metabolomics screening of 1363 metabolites in apparently healthy middle aged (age 45-65) individuals (N = 496) in whom HTGC was measured by proton magnetic resonance spectroscopy. An atlas of metabolite-HTGC associations, based on univariate results, was created using correlation-based Gaussian graphical model (GGM) and genome scale metabolic model network analyses. Pathways associated with the clinical prognosis marker fibrosis 4 (FIB-4) index were tested using a closed global test. RESULTS: Our analyses revealed that 118 metabolites were univariately associated with HTGC (p-value <6.59 × 10-5 ), including 106 endogenous, 1 xenobiotic and 11 partially characterized/uncharacterized metabolites. These associations were mapped to several biological pathways including branched amino acids (BCAA), diglycerols, sphingomyelin, glucosyl-ceramide and lactosyl-ceramide. We also identified a novel possible HTGC-related pathway connecting glutamate, metabolonic lactone sulphate and X-15245 using the GGM network. These pathways were confirmed to be associated with the FIB-4 index as well. The full interactive metabolite-HTGC atlas is provided online: https://tofaquih.github.io/AtlasLiver/. CONCLUSIONS: The combined network and pathway analyses indicated extensive associations between BCAA and the lipids pathways with HTGC and the FIB-4 index. Moreover, we report a novel pathway glutamate-metabolonic lactone sulphate-X-15245 with a potential strong association with HTGC. These findings can aid elucidating HTGC metabolomic profiles and provide insight into novel drug targets for fibrosis-related outcomes.

Alterations in ether lipid metabolism and the consequences for the mouse lipidome.

Alkylglycerol monooxygenase (AGMO) and plasmanylethanolamine desaturase (PEDS1) are enzymes involved in ether lipid metabolism. While AGMO degrades plasmanyl lipids by oxidative cleavage of the ether bond, PEDS1 exclusively synthesizes a specific subclass of ether lipids, the plasmalogens, by introducing a vinyl ether double bond into plasmanylethanolamine phospholipids. Ether lipids are characterized by an ether linkage at the sn-1 position of the glycerol backbone and they are found in membranes of different cell types. Decreased plasmalogen levels have been associated with neurological diseases like Alzheimer's disease. Agmo-deficient mice do not present an obvious phenotype under unchallenged conditions. In contrast, Peds1 knockout mice display a growth phenotype. To investigate the molecular consequences of Agmo and Peds1 deficiency on the mouse lipidome, five tissues from each mouse model were isolated and subjected to high resolution mass spectrometry allowing the characterization of up to 2013 lipid species from 42 lipid subclasses. Agmo knockout mice moderately accumulated plasmanyl and plasmenyl lipid species. Peds1-deficient mice manifested striking changes characterized by a strong reduction of plasmenyl lipids and a concomitant massive accumulation of plasmanyl lipids resulting in increased total ether lipid levels in the analyzed tissues except for the class of phosphatidylethanolamines where total levels remained remarkably constant also in Peds1 knockout mice. The rate-limiting enzyme in ether lipid metabolism, FAR1, was not upregulated in Peds1-deficient mice, indicating that the selective loss of plasmalogens is not sufficient to activate the feedback mechanism observed in total ether lipid deficiency.

Clustered Mendelian randomization analyses identify distinct and opposing pathways in the association between genetically influenced insulin-like growth factor-1 and type 2 diabetes mellitus.

BACKGROUND: There is inconsistent evidence for the causal role of serum insulin-like growth factor-1 (IGF-1) concentration in the pathogenesis of human age-related diseases such as type 2 diabetes (T2D). Here, we investigated the association between IGF-1 and T2D using (clustered) Mendelian randomization (MR) analyses in the UK Biobank. METHODS: We conducted Cox proportional hazard analyses in 451 232 European-ancestry individuals of the UK Biobank (55.3% women, mean age at recruitment 56.6 years), among which 13 247 individuals developed type 2 diabetes during up to 12 years of follow-up. In addition, we conducted two-sample MR analyses based on independent single nucleotide polymorphisms (SNPs) associated with IGF-1. Given the heterogeneity between the MR effect estimates of individual instruments (P-value for Q statistic = 4.03e-145), we also conducted clustered MR analyses. Biological pathway analyses of the identified clusters were performed by over-representation analyses. RESULTS: In the Cox proportional hazard models, with IGF-1 concentrations stratified in quintiles, we observed that participants in the lowest quintile had the highest relative risk of type 2 diabetes [hazard ratio (HR): 1.31; 95% CI: 1.23-1.39). In contrast, in the two-sample MR analyses, higher genetically influenced IGF-1 was associated with a higher risk of type 2 diabetes. Based on the heterogeneous distribution of MR effect estimates of individual instruments, six clusters of genetically determined IGF-1 associated either with a lower or a higher risk of type 2 diabetes were identified. The main clusters in which a higher IGF-1 was associated with a lower risk of type 2 diabetes consisted of instruments mapping to genes in the growth hormone signalling pathway, whereas the main clusters in which a higher IGF-1 was associated with a higher risk of type 2 diabetes consisted of instruments mapping to genes in pathways related to amino acid metabolism and genomic integrity. CONCLUSIONS: The IGF-1-associated SNPs used as genetic instruments in MR analyses showed a heterogeneous distribution of MR effect estimates on the risk of type 2 diabetes. This was likely explained by differences in the underlying molecular pathways that increase IGF-1 concentration and differentially mediate the effects of IGF-1 on type 2 diabetes.

Angiopoietin-like 4 governs diurnal lipoprotein lipase activity in brown adipose tissue.

OBJECTIVE: Brown adipose tissue (BAT) burns fatty acids (FAs) to produce heat, and shows diurnal oscillation in glucose and triglyceride (TG)-derived FA-uptake, peaking around wakening. Here we aimed to gain insight in the diurnal regulation of metabolic BAT activity. METHODS: RNA-sequencing, chromatin immunoprecipitation (ChIP)-sequencing, and lipidomics analyses were performed on BAT samples of wild type C57BL/6J mice collected at 3-hour intervals throughout the day. Knockout and overexpression models were used to study causal relationships in diurnal lipid handling by BAT. RESULTS: We identified pronounced enrichment of oscillating genes involved in extracellular lipolysis in BAT, accompanied by oscillations of FA and monoacylglycerol content. This coincided with peak lipoprotein lipase (Lpl) expression, and was predicted to be driven by peroxisome proliferator-activated receptor gamma (PPARγ) activity. ChIP-sequencing for PPARγ confirmed oscillation in binding of PPARγ to Lpl. Of the known LPL-modulators, angiopoietin-like 4 (Angptl4) showed the largest diurnal amplitude opposite to Lpl, and both Angptl4 knockout and overexpression attenuated oscillations of LPL activity and TG-derived FA-uptake by BAT. CONCLUSIONS: Our findings highlight involvement of PPARγ and a crucial role of ANGPTL4 in mediating the diurnal oscillation of TG-derived FA-uptake by BAT, and imply that time of day is essential when targeting LPL activity in BAT to improve metabolic health.

Deleterious variants in CRLS1 lead to cardiolipin deficiency and cause an autosomal recessive multi-system mitochondrial disease.

Mitochondrial diseases are a group of inherited diseases with highly varied and complex clinical presentations. Here, we report four individuals, including two siblings, affected by a progressive mitochondrial encephalopathy with biallelic variants in the cardiolipin biosynthesis gene CRLS1. Three affected individuals had a similar infantile presentation comprising progressive encephalopathy, bull's eye maculopathy, auditory neuropathy, diabetes insipidus, autonomic instability, cardiac defects and early death. The fourth affected individual presented with chronic encephalopathy with neurodevelopmental regression, congenital nystagmus with decreased vision, sensorineural hearing loss, failure to thrive and acquired microcephaly. Using patient-derived fibroblasts, we characterized cardiolipin synthase 1 (CRLS1) dysfunction that impaired mitochondrial morphology and biogenesis, providing functional evidence that the CRLS1 variants cause mitochondrial disease. Lipid profiling in fibroblasts from two patients further confirmed the functional defect demonstrating reduced cardiolipin levels, altered acyl-chain composition and significantly increased levels of phosphatidylglycerol, the substrate of CRLS1. Proteomic profiling of patient cells and mouse Crls1 knockout cell lines identified both endoplasmic reticular and mitochondrial stress responses, and key features that distinguish between varying degrees of cardiolipin insufficiency. These findings support that deleterious variants in CRLS1 cause an autosomal recessive mitochondrial disease, presenting as a severe encephalopathy with multi-systemic involvement. Furthermore, we identify key signatures in cardiolipin and proteome profiles across various degrees of cardiolipin loss, facilitating the use of omics technologies to guide future diagnosis of mitochondrial diseases.

Adverse generational changes in obesity development converge at midlife without increased cardiometabolic risk.

OBJECTIVE: Obesity is becoming a global public health problem, but it is unclear how it impacts different generations over the life course. Here, a descriptive analysis of the age-related changes in anthropometric measures and related cardiometabolic risk factors across different generations was performed. METHODS: The development of anthropometric measures and related cardiometabolic risk factors was studied during 26 years of follow-up in the Doetinchem Cohort Study (N = 6,314 at baseline). All analyses were stratified by sex and generation, i.e., 10-year age groups (20-29, 30-39, 40-49, and 50-59 years) at baseline. Generalized estimating equations were used to test for generational differences. RESULTS: Weight, BMI, waist circumference, and prevalence of overweight and obesity were higher, in general, in the younger generations during the first 10 to 15 years of follow-up. From age 50 to 59 years onward, these measures converged in all generations of men and women. Among cardiometabolic risk factors, only type 2 diabetes showed an unfavorable shift between the two oldest generations of men. CONCLUSIONS: It was observed that, compared with the older generations, the younger generations had obesity at an earlier age but did not reach higher levels at midlife and beyond. This increased exposure to obesity was not (yet) associated with increased prevalence of cardiometabolic risk factors.

Genome-Wide Association Study of the Postprandial Triglyceride Response Yields Common Genetic Variation in LIPC (Hepatic Lipase).

BACKGROUND: The increase in serum triglyceride (TG) concentrations in response to a meal is considered a risk factor for cardiovascular disease. We aimed to elucidate the genetics of the postprandial TG response through genome-wide association studies (GWAS). METHODS: Participants of the NEO (Netherlands Epidemiology of Obesity) study (n=5630) consumed a liquid mixed meal after an overnight fast. GWAS of fasting and postprandial serum TG at 150 minutes were performed. To identify genetic variation of postprandial TG independent of fasting TG, we calculated the TG response at 150 minutes by the residuals of a nonlinear regression that predicted TG at 150 minutes as a function of fasting TG. Association analyses were adjusted for age, sex, and principal components in a linear regression model. Next, using the identified variants as determinants, we performed linear regression analyses on the residuals of the postprandial response of 149 nuclear magnetic resonance-based metabolite measures. RESULTS: GWAS of fasting TG and postprandial serum TG at 150 minutes resulted in completely overlapping loci, replicating previous GWAS. From GWAS of the TG response, we identified rs7350789-A (allele frequency=0.36), mapping to hepatic lipase (LIPC), to be associated with a smaller increase in TG concentrations at 150 minutes (ß=-0.11; P-value=5.1×10-8). Rs7350789-A was associated with responses of 33 metabolite measures (P-value <1.34×10-3), mainly smaller increases of the TG-component in almost all HDL (high-density lipoprotein) subparticles (HDL-TG), a smaller decrease of HDL diameter and smaller increases of most components of VLDL (very low density lipoprotein) subparticles. CONCLUSIONS: GWAS of the TG response identified a variant near LIPC as a main contributor to postprandial TG metabolism independent of fasting TG concentrations, resulting in smaller increases of HDL-TG and VLDL subparticles.

The effect of mirabegron on energy expenditure and brown adipose tissue in healthy lean South Asian and Europid men.

AIM: To compare the effects of cold exposure and the ß3-adrenergic receptor agonist mirabegron on plasma lipids, energy expenditure and brown adipose tissue (BAT) activity in South Asians versus Europids. MATERIALS AND METHODS: Ten lean Dutch South Asian (aged 18-30 years; body mass index [BMI] 18-25 kg/m2 ) and 10 age- and BMI-matched Europid men participated in a randomized, double-blinded, cross-over study consisting of three interventions: short-term (~ 2 hours) cold exposure, mirabegron (200 mg one dose p.o.) and placebo. Before and after each intervention, we performed lipidomic analysis in serum, assessed resting energy expenditure (REE) and skin temperature, and measured BAT fat fraction by magnetic resonance imaging. RESULTS: In both ethnicities, cold exposure increased the levels of several serum lipid species, whereas mirabegron only increased free fatty acids. Cold exposure increased lipid oxidation in both ethnicities, while mirabegron increased lipid oxidation in Europids only. Cold exposure and mirabegron enhanced supraclavicular skin temperature in both ethnicities. Cold exposure decreased BAT fat fraction in both ethnicities. After the combination of data from both ethnicities, mirabegron decreased BAT fat fraction compared with placebo. CONCLUSIONS: In South Asians and Europids, cold exposure and mirabegron induced beneficial metabolic effects. When combining both ethnicities, cold exposure and mirabegron increased REE and lipid oxidation, coinciding with a higher supraclavicular skin temperature and lower BAT fat fraction.

Integration of epidemiologic, pharmacologic, genetic and gut microbiome data in a drug-metabolite atlas.

Progress in high-throughput metabolic profiling provides unprecedented opportunities to obtain insights into the effects of drugs on human metabolism. The Biobanking BioMolecular Research Infrastructure of the Netherlands has constructed an atlas of drug-metabolite associations for 87 commonly prescribed drugs and 150 clinically relevant plasma-based metabolites assessed by proton nuclear magnetic resonance. The atlas includes a meta-analysis of ten cohorts (18,873 persons) and uncovers 1,071 drug-metabolite associations after evaluation of confounders including co-treatment. We show that the effect estimates of statins on metabolites from the cross-sectional study are comparable to those from intervention and genetic observational studies. Further data integration links proton pump inhibitors to circulating metabolites, liver function, hepatic steatosis and the gut microbiome. Our atlas provides a tool for targeted experimental pharmaceutical research and clinical trials to improve drug efficacy, safety and repurposing. We provide a web-based resource for visualization of the atlas (http://bbmri.researchlumc.nl/atlas/).

A metabolic profile of all-cause mortality risk identified in an observational study of 44,168 individuals.

Predicting longer-term mortality risk requires collection of clinical data, which is often cumbersome. Therefore, we use a well-standardized metabolomics platform to identify metabolic predictors of long-term mortality in the circulation of 44,168 individuals (age at baseline 18-109), of whom 5512 died during follow-up. We apply a stepwise (forward-backward) procedure based on meta-analysis results and identify 14 circulating biomarkers independently associating with all-cause mortality. Overall, these associations are similar in men and women and across different age strata. We subsequently show that the prediction accuracy of 5- and 10-year mortality based on a model containing the identified biomarkers and sex (C-statistic = 0.837 and 0.830, respectively) is better than that of a model containing conventional risk factors for mortality (C-statistic = 0.772 and 0.790, respectively). The use of the identified metabolic profile as a predictor of mortality or surrogate endpoint in clinical studies needs further investigation.

Reduced Feed Intake, Rather than Increased Energy Losses, Explains Variation in Growth Rates of Normal-Birth-Weight Piglets.

Background: Substantial variation in growth rates exists in normal-birth-weight piglets, possibly due to differences in energy efficiency. Within this population, slow growth rates are associated with reduced insulin sensitivity. Slowly digestible starch (SDS) may improve growth efficiency in slowly growing pigs, because it reduces postprandial blood glucose. Objective: The aim of this study was to investigate maintenance energy requirements and efficiency of energy used for growth (incremental energy efficiency) of slow-growing or fast-growing piglets (SG-pigs and FG-pigs, respectively) with equal birth weight that were fed either an SDS or a rapidly digestible-starch (RDS) diet. Methods: Sixteen groups of either five 10-wk-old SG-pigs (mean ± SD: 11.3 ± 1.4 kg) or FG-pigs (15.1 ± 1.7 kg) were housed in climate respiration chambers and fed diets containing 40% RDS or SDS for 2 wk. In week 1, feed was available ad libitum. In week 2, feed supply was restricted to 65% of the observed weekly averaged feed intake [kJ · kg body weight (BW)-0.6 · d-1] in week 1. After week 2, pigs were feed deprived for 24 h, after which heat production was determined. Energy balances, apparent total tract digestibility (ATTD), and incremental energy efficiencies were calculated and analyzed using a general linear model. Results: Gross energy intake (kJ · kg BW-0.6 · d-1) was 4% greater (P = 0.047) for FG-pigs than for SG-pigs. ATTD of fat was 6%-units greater (P = 0.003) for RDS-fed than for SDS-fed pigs. Fasting heat production and incremental energy efficiencies did not differ between pig types or diets. Incremental use of metabolizable energy for fat retention was 2% units (P = 0.054) greater for RDS-fed than SDS-fed pigs. Conclusions: A lower energy intake rather than greater maintenance requirements or lower energy efficiency explains the slow growth of SG-pigs. Incremental RDS intake increased fat deposition more than SDS, whereas energy efficiency was not affected. Thus, feeding SDS instead of RDS does not improve growth efficiency but may result in slightly leaner pigs.

Pyruvate dehydrogenase complex plays a central role in brown adipocyte energy expenditure and fuel utilization during short-term beta-adrenergic activation.

Activation of brown adipose tissue (BAT) contributes to total body energy expenditure through energy dissipation as heat. Activated BAT increases the clearance of lipids and glucose from the circulation, but how BAT accommodates large influx of multiple substrates is not well defined. The purpose of this work was to assess the metabolic fluxes in brown adipocytes during ß3-adrenergic receptor (ß3-AR) activation.T37i murine preadipocytes were differentiated into brown adipocytes and we used Seahorse respirometry employing a set of specific substrate inhibitors in the presence or absence of ß3-AR agonist CL316,243. The main substrate used by these brown adipocytes were fatty acids, which were oxidized equally during activation as well as during resting condition. [U-13C]-glucose tracer-based metabolomics revealed that the flux through the TCA cycle was enhanced and regulated by pyruvate dehydrogenase (PDH) activity. Based on 13C-tracer incorporation in lipids, it appeared that most glucose was oxidized via TCA cycle activity, while some was utilized for glycerol-3-phosphate synthesis to replenish the triglyceride pool. Collectively, we show that while fatty acids are the main substrates for oxidation, glucose is also oxidized to meet the increased energy demand during short term ß3-AR activation. PDH plays an important role in directing glucose carbons towards oxidation.

Postprandial metabolite profiles associated with type 2 diabetes clearly stratify individuals with impaired fasting glucose.

INTRODUCTION: Fasting metabolite profiles have been shown to distinguish type 2 diabetes (T2D) patients from normal glucose tolerance (NGT) individuals. OBJECTIVES: We investigated whether, besides fasting metabolite profiles, postprandial metabolite profiles associated with T2D can stratify individuals with impaired fasting glucose (IFG) by their similarities to T2D. METHODS: Three groups of individuals (age 45-65 years) without any history of IFG or T2D were selected from the Netherlands Epidemiology of Obesity study and stratified by baseline fasting glucose concentrations (NGT (n = 176), IFG (n = 186), T2D (n = 171)). 163 metabolites were measured under fasting and postprandial states (150 min after a meal challenge). Metabolite profiles specific for a high risk of T2D were identified by LASSO regression for fasting and postprandial states. The selected profiles were utilised to stratify IFG group into high (T2D probability ≥ 0.7) and low (T2D probability ≤ 0.5) risk subgroups. The stratification performances were compared with clinically relevant metabolic traits. RESULTS: Two metabolite profiles specific for T2D (nfasting = 12 metabolites, npostprandial = 4 metabolites) were identified, with all four postprandial metabolites also being identified in the fasting state. Stratified by the postprandial profile, the high-risk subgroup of IFG individuals (n = 72) showed similar glucose concentrations to the low-risk subgroup (n = 57), yet a higher BMI (difference: 3.3 kg/m2 (95% CI 1.7-5.0)) and postprandial insulin concentrations (21.5 mU/L (95% CI 1.8-41.2)). CONCLUSION: Postprandial metabolites identified T2D patients as good as fasting metabolites and exhibited enhanced signals for IFG stratification, which offers a proof of concept that metabolomics research should not focus on the fasting state alone.

A Mendelian Randomization Study of Metabolite Profiles, Fasting Glucose, and Type 2 Diabetes.

Mendelian randomization (MR) provides us the opportunity to investigate the causal paths of metabolites in type 2 diabetes and glucose homeostasis. We developed and tested an MR approach based on genetic risk scoring for plasma metabolite levels, utilizing a pathway-based sensitivity analysis to control for nonspecific effects. We focused on 124 circulating metabolites that correlate with fasting glucose in the Erasmus Rucphen Family (ERF) study (n = 2,564) and tested the possible causal effect of each metabolite with glucose and type 2 diabetes and vice versa. We detected 14 paths with potential causal effects by MR, following pathway-based sensitivity analysis. Our results suggest that elevated plasma triglycerides might be partially responsible for increased glucose levels and type 2 diabetes risk, which is consistent with previous reports. Additionally, elevated HDL components, i.e., small HDL triglycerides, might have a causal role of elevating glucose levels. In contrast, large (L) and extra large (XL) HDL lipid components, i.e., XL-HDL cholesterol, XL-HDL-free cholesterol, XL-HDL phospholipids, L-HDL cholesterol, and L-HDL-free cholesterol, as well as HDL cholesterol seem to be protective against increasing fasting glucose but not against type 2 diabetes. Finally, we demonstrate that genetic predisposition to type 2 diabetes associates with increased levels of alanine and decreased levels of phosphatidylcholine alkyl-acyl C42:5 and phosphatidylcholine alkyl-acyl C44:4. Our MR results provide novel insight into promising causal paths to and from glucose and type 2 diabetes and underline the value of additional information from high-resolution metabolomics over classic biochemistry.

Weight loss induced by very low calorie diet is associated with a more beneficial systemic inflammatory profile than by Roux-en-Y gastric bypass.

INTRODUCTION: Weight loss interventions such as Roux-en-Y gastric bypass (RYGB) and very low calorie diets (VLCD) lead to improvement of glucose metabolism in obese individuals with type-2 diabetes. Weight loss can also positively influence the unfavorable inflammatory profile associated with obesity. However, a direct comparison of the effect of VLCD and RYGB on systemic inflammation is lacking. METHODS: Systemic inflammation was investigated in age- and BMI-matched morbidly obese T2DM women by determining the number and activation- or memory status of peripheral blood leukocytes by flow cytometry, in addition to measuring circulating levels of cytokines and CRP. Systemic inflammation was assessed one month before and three months after RYGB (n=15) or VLCD (n=12). An age matched group of lean women (n=12) was studied as control group. RESULTS: Three months after the intervention, CRP and leptin levels were reduced whereas adiponectin levels were increased both by RYGB and VLCD. TNF-α levels were increased by RYGB, but reduced by VLCD. IL-2 and IL-6 levels were reduced and IL-4 levels were increased by VLCD but not affected by RYGB. The number of activated peripheral cytotoxic T (CD8+CD25+) and B (CD19+CD38+) cells was significantly higher after RYGB than after VLCD. CONCLUSION: In conclusion, RYGB and VLCD have differential effects on the activation status of peripheral leukocytes and levels of cytokines in obese women with T2DM, despite comparable weight loss three months after the intervention. VLCD seems to have more favorable effects on the inflammatory profile as compared to RYGB.

Reanalysis of mGWAS results and in vitro validation show that lactate dehydrogenase interacts with branched-chain amino acid metabolism.

The assignment of causative genes to noncoding variants identified in genome-wide association studies (GWASs) is challenging. We show how combination of knowledge from gene and pathway databases and chromatin interaction data leads to reinterpretation of published quantitative trait loci for blood metabolites. We describe a previously unidentified link between the rs2403254 locus, which is associated with the ratio of 3-methyl-2-oxobutanoate and alpha-hydroxyisovalerate levels, and the distal LDHA gene. We confirmed that lactate dehydrogenase can catalyze the conversion between these metabolites in vitro, suggesting that it has a role in branched-chain amino acid metabolism. Examining datasets from the ENCODE project we found evidence that the locus and LDHA promoter physically interact, showing that LDHA expression is likely under control of distal regulatory elements. Importantly, this discovery demonstrates that bioinformatic workflows for data integration can have a vital role in the interpretation of GWAS results.

FluxModeCalculator: an efficient tool for large-scale flux mode computation.

UNLABELLED: Elementary flux mode (EFM) analysis is a powerful technique for determining the metabolic capacities and robustness of stoichiometric networks. Recently, several improvements have been made to the algorithm for enumerating the EFMs, making the study of large models possible. However, currently these tools require high performance workstations to perform large-scale EFM computations, thus limiting their applicability. We developed a more time and memory efficient implementation of the algorithm for EFM enumeration in MATLAB, called FluxModeCalculator, which enables large-scale EFM computation on ordinary desktop computers. AVAILABILITY AND IMPLEMENTATION: FluxModeCalculator is open source and freely available under the terms of the GNU General Public License v3.0 at http://www.lumc.nl/jan-bert-van-klinken CONTACT: j.b.van_klinken@lumc.nl SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genome-wide association study identifies novel genetic variants contributing to variation in blood metabolite levels.

Metabolites are small molecules involved in cellular metabolism, which can be detected in biological samples using metabolomic techniques. Here we present the results of genome-wide association and meta-analyses for variation in the blood serum levels of 129 metabolites as measured by the Biocrates metabolomic platform. In a discovery sample of 7,478 individuals of European descent, we find 4,068 genome- and metabolome-wide significant (Z-test, P < 1.09 × 10(-9)) associations between single-nucleotide polymorphisms (SNPs) and metabolites, involving 59 independent SNPs and 85 metabolites. Five of the fifty-nine independent SNPs are new for serum metabolite levels, and were followed-up for replication in an independent sample (N = 1,182). The novel SNPs are located in or near genes encoding metabolite transporter proteins or enzymes (SLC22A16, ARG1, AGPS and ACSL1) that have demonstrated biomedical or pharmaceutical importance. The further characterization of genetic influences on metabolic phenotypes is important for progress in biological and medical research.